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ABSTRACT 

We analyze several recently detected gamma-ray bursts (GRBs) with late X-ray flares in the context 
of late internal shock and late external shock models. We And that the X-ray flares in GRB 050421 
and GRB 050502B originate from late internal shocks, while the main X-ray flares in GRB 050406 and 
GRB 050607 may arise from late external shocks. Under the assumption that the central engine has 
two periods of activities, we get four basic types of X-ray light curves. The classification of these types 
depends on which period of activities produces the prompt gamma-ray emission (Type 1 and Type 2: 
the earlier period; Type 3 and Type 4: the late period), and on whether the late ejecta catching up with 
the early ejecta happens earlier than the deceleration of the early ejecta (Type 1 and Type 3) or not 
(Type 2 and Type 4). We And that the X-ray flare caused by a late external shock is a special case of 
Type 1. Our analysis reveals that the X-ray light curves of GRBs 050406, 050421, and 050607 can be 
classified as Type 1, while the X-ray light curve of GRB 050502B is classified as Type 2. However, the 
X-ray light curve of GRB 050406 is also likely to be Type 2. We also predict a long-lag short-lived X-ray 
flare caused by the inner external shock, which forms when a low baryon-loading long-lag late ejecta 
decelerates in the non-relativistic tail of an outer external shock driven by an early ejecta. 

Subject headings: gamma rays: bursts—X-rays: flares— shocks: relativistic 


1. INTRODUCTION 

In the pre-Swift era, only a few early optical afterglows 
of gamma-ray bursts (GRBs) were detected, while the ob¬ 
servations of X-ray afterglows usually started at ~ 10"^ sec¬ 
onds after the prompt trigger. Since early afterglows con¬ 
tain important information about GRB central engines, 
understanding the early afterglows is therefore one of the 
most interesting scientific goals of the NASA’s Swift satel¬ 
lite (Gehrels et al. 2004). After about one year of oper¬ 
ations, a significant fraction of well localized GRBs have 
not been detected with early optical emissions down to 
moderate limiting magnitudes, although the UV/Optical 
Telescope (UVOT) on board Swift slewed to the error cir¬ 
cle of position quickly after the burst (Roming et al. 2005). 
On the other hand, observations by the X-ray Telescope 
(XRT) on board Swift have revealed several new features of 
X-ray emissions. First, steep declines in some X-ray light 
curves at the transition from prompt phase to afterglow 
phase have been discovered, which are interpreted as tail 
emissions of prompt GRBs (Tagliaferri et al. 2005; Zhang 
et al. 2005; Nousek et al. 2005). Secondly, a portion 
of early X-ray afterglows have shallow-than-normal tem¬ 
poral decays before they enter “normal” decaying phase 
(Nousek et al. 2005; Zhang et al. 2005). Thirdly, late X- 
ray flares have been observed in several GRBs (Burrows 
et al. 2005a; see also Burrows et al. 2005b for a review). 
Piro et al. (2005) also reported that X-ray flares were dis¬ 
covered in a few GRBs (e.g., GRBs 011121 and 011211) 
by the Italian-Dutch BeppoSAX satellite. GRBs are usu¬ 
ally divided into two main classes by their durations and 
spectral hardness ratios, i.e., long GRBs have durations 


larger than 2 seconds and softer spectra while short GRBs 
have durations shorter than 2 seconds and harder spectra 
(Kouveliotou et al. 1993). Up to now, both long GRBs, 
including X-ray rich bursts and a high redshift burst (GRB 
050904), and short GRBs have been found with late X-ray 
flares (Galli & Piro 2005; Watson et al. 2005; Fox et al. 
2005; Barthelmy et al. 2005). 

One leading explanation for the early shallow decaying 
X-ray afterglows is that the forward external shock has 
continuous energy injection from the long active central 
engine (Zhang et al. 2005; Nousek et al. 2005; Dai & Lu 
1998a; Zhang & Meszaros 2002; Dai 2004; Wei, Yan & Fan 
2005). Another leading explanation is that the Lorentz 
factor of the GRB ejecta has a distribution shaped by the 
central engine so that the behind slower material catches 
up with the ahead faster material when the latter is decel¬ 
erated in the circum-burst medium, acting as a continuous 
energy injection (Zhang et al. 2005; Nousek et al. 2005; 
Rees & Meszaros 1998; Sari & Meszaros 2000; Zhang & 
Meszaros 2002; Granot & Kumar 2005). In this explana¬ 
tion the central engine does not need to be active for a long 
time (Granot & Kumar 2005). It is also possible that the 
early shallow decaying X-ray afterglows are caused by late 
continuous energy releases if the GRB ejecta is initially 
dominated by Poynting flux (Zhang & Kobayashi 2004). 

As for X-ray flares, both late external shock model and 
late internal shock model have been proposed (Piro et al. 
2005; Burrows et al. 2005; Fan & Wei 2005; Zhang et al. 
2005). Piro et al. (2005) and Galli & Piro (2005) argued 
for the late external shock model because they found there 
is no obvious spectral evolution in some X-ray flares, while 
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Burrows et al. (2005b) preferred the late internal shock 
model based on several other X-ray flares which showed 
opposite features. These two conclusions are only based 
on their qualitative and empirical analysis. In this paper 
we quantitatively compare the X-ray flares with the pre¬ 
dictions of the above two models. In Section 2 we describe 
the intrinsic X-ray light curves from the internal and exter¬ 
nal shocks. Then we analyze four GRBs with X-ray flares 
in Section 3 by comparing with the “delayed” (relative to 
the GRB trigger) intrinsic light curves of the late internal 
shock and late external shock. In Section 4 we establish 
four basic types of X-ray light curves, assuming that the 
central engine has two periods of activities. We also clas¬ 
sify the four GRBs by these types. Our discussion and 
conclusions about X-ray flares are presented in Section 5. 

2. INTRINSIC LIGHT CURVES OF INTERNAL AND 
EXTERNAL SHOCK EMISSIONS 

Suppose that an intrinsic light curve of the inter¬ 
nal/external shock emissions comprises an initial rise and 
a subsequent decay, which can approximately be described 
by the following broken power-law function of the ob¬ 
server’s time, 

p _ \ Pv,pk{t/tbT'^, t<tb, 

" 1 i^ppfe(fAh)“G t>tb, 

where oi (> 0), a 2 (< 0) are the temporal indices before 
and after the peak time tb, and F^^pk is the peak flux den¬ 
sity correspondingly. The FWHM width of the rising and 
falling time scale of the light curve are Str = {l — 2~^'°‘^)tb 
and 6tf = (2“^/“^ — 1)4, respectively. 

In the internal shock model, two shells with a difference 
in Lorentz factor (T) AT ~ T ^ 1 and a time delay At 
by the central engine collide at the radius i?int ~ 2r^cAt 
(Paczynski & Xu 1994). During the interaction of these 
two shells, the observed emission flux rises rapidly. The 
high time-resolution observations of single-pulse GRBs 
(e.g. GRB 971208, Connaughton et al. 1997) and main 
long pulses in usual GRBs reveal that the emission rises 
linearly with time, i.e., ai ~ 1. One should keep in mind 
that this value is just empirical. The falling behavior of 
the light curve is mainly attributed to high-latitude emis¬ 
sion (also known as tail emission) of the fireball, if the 
width of the fireball can be neglected (Fenimore, Madras 
& Nayakshin 1996; Kumar & Panaitescu 2000). Let the 
co-moving flux density F)(, oc be uniform in the fire¬ 
ball, the observed flux density is F^, oc Fl^,'D'^{dT,/dt), 
where T> = l/r(l — cosd) is the Doppler factor and 
the surface area S satisfies dS = sin 9d9 oc d9^ 

(Panaitescu et al. 2005). The time delay between pho¬ 
tons emitted at an angle 9 and those at the line of sight 
(LOS, 9 = 0)is t — tb = i?int(l — cosd)/c « j^c = 
(dF)^At, which further leads to dYi/dt =constant and 
V « 2r/[(t — tb)/At -I- 1]. Since F = v/V, the theo¬ 
retical falling behavior of the internal shock emission is 
F^ oc oc [(t —4)/At-1-1]“^^“^) (Kumar & Panaitescu 

2000), or exactly 

F, = F,^pk[{t-tb)IAt + l]-^'^-P\ tb<t. ( 2 ) 

The initial {tb < t < At + tb) decaying part of the light 
curve can be approximated by an exponential function of 


time, i.e., F^ oc , where r = At/{2 — /3) ~ At/3. 
This behavior has been observed in most of GRBs. At 
late times (t At -I- tf,), the tail emission decays as a 
power law function of time, Fi, oc which has been 

confirmed by recent X-ray observations conducted dur¬ 
ing the prompt GRB phase by the Swift satellite (Tagli- 
aferri et al. 2005; Nousek et al. 2005). The decaying 
time scale (FWHM) of internal shock emission is therefore 
dtf = — l]At ~ (ln2)T, not directly related to tb- 

It should be noted that the observed Str/6tf is ~ 0.3 — 0.5 
in prompt GRBs (Norris et al. 1996). The temporal index 
02 = d\nF„/d\nt = (/3 — 2)t/{t — 4 + At), ranges from 
— (2 — j3)tb/At at t ~ 4 to —(2 — j3) when t^ tb + At. In 
general, the light curve of internal shock emission can be 
roughly depicted by eq. (1), especially for the rising and 
very late decaying portions with ai ~ 1 and a 2 = -(2-/3). 

In the external shock model, if the radial width of the 
GRB ejecta is smaller than a critical value (thin shell 
case), the peak time 4 is regarded as the deceleration time 
4ec of the external shock sweeping into the circum-burst 
medium, driven by the ejecta. On the other hand, if the 
width of the ejecta is larger than that critical value (thick 
shell case), 4 is the observed crossing time of the reverse 
shock through the original GRB ejecta, which is about the 
width divided by the speed of light (Sari & Piran 1999). 
Since the typical Lorentz factor of electrons in the forward 
shock is much larger than that in the reverse shock, the 
X-ray emission in most common cases is dominated by the 
forward shock (Sari & Piran 1999; see, however, Fan & 
Wei 2005 and Kobayashi et al. 2005, for opposite cases 
under some extreme assumptions). 

If the circum-burst environment is an interstellar 
medium (ISM), the reverse shock velocity is Newtonian 
or relativistic depending on whether the initial Lorentz 
factor Fi is smaller or larger than the critical value Fc « 
200if53®nQ where n = no cm“^ is the ISM den¬ 

sity, E = 10^^i?53 erg and A^ = 10^^Aip2 cm are the 
isotropic-equivalent kinetic energy and initial width of the 
ejecta (Zhang, Kobayashi & Meszaros 2003). For the New¬ 
tonian reverse shock (NRS) case, the Lorentz factor of 
the forward shock roughly equals its initial value before 
the reverse shock crosses the ejecta, F « F^. The mini¬ 
mum Lorentz factor of shock-accelerated electrons and co¬ 
moving magnetic field in the forward shock scale as ■jm oc F 
and B oc Fn^/^ respectively, while the cooling Lorentz fac¬ 
tor of electrons in the dynamical time scale is jc oc l/tT^n 
(Sari, Piran & Narayan 1998). The characteristic fre¬ 
quencies and maximum flux density of the synchrotron 
radiation from the forward shock are oc Fy^i? ex 
Vc cx Fy/R oc and Fly, max oc nR^TB oc t^ (Sari & Pi¬ 
ran 1999). For the relativistic reverse shock (RRS) case, 
the radius and Lorentz factor of the forward shock evolve 
as i? oc and F oc t~^l'^, respectively (Kobayashi 2000). 
The characteristic frequencies and peak flux density in 
this case are v^a oc t~^, Vc oc t~^, and F,y,max oc t. Af¬ 
ter the reverse shock crosses the ejecta, the forward shock 
emission is characterized by Vm oc Vc oc and 

Fi/,max oc (Sari, Piran & Narayan 1998). Table 1 sum¬ 
marizes the temporal indices ai and 02 , and the corre¬ 
sponding spectral index /3 (F^, oc v^) oi the forward shock 
emission in ISM. We can see that the flux density rises 
rapidly in the NRS case, while it rises quite slowly in the 
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RRS case. Note that in the RRS case if the crossing times 
tm (when Vm = smd tc (when v,, = v) are both earlier 
than the time when the reverse shock crosses the ejecta, 
the peak time th should be the maximum of tm and tc- 
In such a special case, ai = 1/2 {tc < tm) or (3 — p)/2 
{tc > tm), 02 equals to —{p — 2)f2 initially and then to 
— (3p —2)/4. The ratio of the rising to falling times 6tr/Stf 
in this case ranges from ~ 10“^ to ~ 1, and the spectral 
index changes from —1/2 or —{p— l)/2 to —p/2 around 
tb- 

If the circum-burst environment is a stellar wind with 
density n = AR~^ (A = 3 x 10^®A* cm“^ is the wind 
parameter, Chevalier & Li 2000; Dai & Lu 1998b), the re¬ 
verse shock velocity is Newtonian or relativistic depending 
on whether the initial Lorentz factor T i is smaller or larger 
than the critical value Tc « (Zou, Wu 

& Dai 2005). However, the Lorentz factor L of the forward 
shock remains to be a constant before the reverse shock 
crosses the ejecta, although L « for the NRS case while 
r ^ Li for the RRS case (Chevalier & Li 2000; Wu et 
al. 2003; Zou, Wu & Dai 2005). Since R (x t, n cc t~'^, 
B (X 7m oc and 7 c oc t, we obtain Vm oc t~^, 
Vc (X t and Tly,max oc After the reverse shock crosses 
the ejecta, the forward shock emission is characterized by 
Vm oc Vc oc and Aly,max oc (Chevalier & Li 

2000). Table 2 summarizes the temporal and spectral in¬ 
dices ai, 02 , and [3 of the forward shock emission in a stel¬ 
lar wind. To reproduce an initial increasing flux density 
in the wind environment one must require Vc < v < I’m- 
If the time when Um equals to v is earlier than the reverse 
shock crosses the ejecta, the peak time tb should be tm- In 
this case, oi = 1/2, a 2 is equal to —{p — 2)/2 initially and 
then to —{3p — 2)/4, Str/Stf ~ 10“^ — 1, and the spectral 
index changes from —1/2 to —p/2 at around tb- 

From the above analysis we can see that the simple in¬ 
ternal or external shock models can not account for the 
observed X-ray flares with fast rising {ai ~ 10) and rapid 
decay (02 ~ —10) behavior.^ We study the late internal 
shock (LIS) model and late external shock (LES) model 
for X-ray flares in detail and make a comparison between 
them case by case in particular flares in the next section. 

3. TIME ZERO POINT EFFECT ON THE DELAYED 
INTERNAL/EXTERNAL SHOCK EMISSIONS 

Assume there exists a time delay to between the trigger 
of a GRB and the beginning of the above light curve as de¬ 
scribed by eq. (1). The influence of this time delay on the 
light curve is known as the time zero point effect (Huang, 
Dai & Lu 2002; Zhang et al. 2005). The physical reason 
and the conditions for the delay within the internal and 
external shock models will be presented in the next sec¬ 
tion. Here we make a quantitative analysis of this effect. 
Due to the time delay, the observed light curve plotted 
with the time zero point chosen as the GRB trigger is 

( Bjcp}^{{t to^/tb] to ^ t <t tpeakj 

Fi, = < F^^pk[{t - to)/4]“L t > tpeak (LES), 

[ F„^pk[{t - tpeak)/At + f (LIS), 

(3) 


where the peak time is changed to tpeak = to + tb- The 
observed temporal index of any segment of the light curve 
becomes (e.g., Huang, Dai & Lu 2002) 


Q^obs — 


dlnFb 

dint 


{P - 2)t 

t - tpeak + At ’ 

at 

t-to’ 


t ^ tpeak (LIS), 
otherwise. 


(4) 

Therefore, the observed temporal indices just before and 
after tpeak are Oi^obs — Clitpeak/^bi CK2,obs — Cl2tpeak/G 
(LES), and a2,obs = {P - 2)tpeak/At (LIS) (see Eig. 1). 
When t tpeak, the to-shifted light curve approaches 
to the intrinsic one. The ratio dt//tpeak, which equals 
_ l)tf,/tpeak in the LES model and (2^A2-/3) _ 
I) At/tpeak in the LIS model, is ~ 0.1 in all of the observed 
X-ray flares (Burrows et al. 2005a, b). This means that 
the time delay to is much larger than tb or At. A conse¬ 
quence of the time zero point effect is that there exists an 
anti-correlation between the two observed quantities aobs 


and dt/tpeak, 


^l,obs — Bi{Strltpeak) , 0^2,obs — B2{Stf /tpeak) , (5) 


where Bi = ai(I —2“^/“^) ranges from ~ 0.4 for ai = 1/2 
to ~ 0.6 for di = 3 and has an upper limit Bi < 0.69, while 
i ?2 = 02 ( 2 “^^“^ — 1) varies from —1 for a 2 = — 1 (typical 
of LES) to —0.8 for a 2 = —3 (typical of LIS, in which case 
02 in eq. [5] can be regarded as P — 2) and has an upper 
limit i ?2 < —0.69. As can be seen in Eig. 1, the ratio of 
the peak flux density directly measured at tpeak to the one 
which is extrapolated from late time {t ^ tpeak) back to 
tpeak, or namely the flare increasing factor, is 


[( At ^ 

v/5-2 //3-2N/5-2 

1 ^ ^peak ^ 

^ Ct2,obs ^ 

\ttb^ 

/ a2 

^ ^ ^peak ^ 

^ 0^2,obs ^ 


LIS, 

LES. 


( 6 ) 


This value is also roughly equal to the ratio of to 
Fi,(2tpeak), as long as to ^ tb or At. Note that for a 
given a 2 .obs, Am has a maximum Am,max « l, 44 -“ 2 ,ob= 
when a 2 ~ 0.37a2,obs- In the late internal shock model, 
adopting /3 — 2 ~ —3.0 and Stf/tpeak ~ 0.1, the flare in¬ 
creasing factor is typically Am ~ 20. Because the late time 
emission is dominated by the external shock, the measured 
Am,obs can be either larger or smaller than this theoret¬ 
ical value, depending on whether the backward extrapo¬ 
lated flux density of the external shock emission at tpeak 
is smaller or larger than that of the LIS emission at this 
time. To reduce the deviation from the theoretical Am 
influenced by the external shock emission, the observed 
Am,obs niay be better determined by the ratio of Eb(tpeak) 
to Eb(2tpeak), rather than the value extrapolated from the 
late time light curve. In the late external shock model, 
however, the measured Am reflects the true ratio. The 
value of 02 can be directly determined by the post-flare 
light curve. Taken typical 02 ~ —1 and Stf/tpeak ~ 0.1, 
we obtain Am ~ 10 for typical LES X-ray flares. 

Hence, once we know ai, 02 , Str, Stf, tpeak and Am,ohs 
of a particular X-ray flare, the intrinsic values of ai, 02 , 


^We do not consider the effect of magnetization of GRB ejecta on the external shock emission, which may change the temporal indices 
moderately as compared to those listed in Tables 1 and 2 (Zhang & Kobayashi 2005). Nevertheless, this effect cannot explain the observed 
rapid rising and decaying features of X-ray flares. 
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tb and to are over-determined. Together with the spectral 
information during the flare, we can tightly constrain the 
X-ray flare models. As the X-ray spectral information is 
often contained in the ratio of the count rate in hard band 
to that in soft band (the so called hardness ratio H/S), 
we convert this ratio to the spectral index (3 with Fig. 2. 
Below we present our case studies for X-ray flares using 
available data in the literature. 

GRB O5O4O6 — The observed quantities are Oi^obs = 
-a2,obs = 6.8, tpeak = 213 S, 6 tr « 5tf = 0.20to;J^tpeak, 
and Am,ohs ~ 6 (Romano et al. 2005; Burrows et al. 2005). 
Using eq (5) we obtain 02 = ~0.56lg;23, consistent with 
0 C 2 ~ —0.88 estimated by eq (6). However, it is impos¬ 
sible for the flare to have a highly symmetric light curve 
with 6tr = Stf and ai,obs = —a2,obs, unless ai = —02 
approaches to infinity. If we do not take seriously on the 
property of the rising behavior and adopt a 2 ~ —0.88, 
then we get 4 = 28 s and to = 185 s. The hardness ra¬ 
tio of this X-ray afterglow evolves from 0.4 ± 0.3 initially 
to 1.3 ± 0.3 at the peak of the hard band (1 — 10 keV) 
light curve, and Anally decreases to a level of ~ 0.7 af¬ 
ter the peak time of the soft band (0.1 — 10 keV) light 
curve. This corresponds to the spectral index (3 being 
—0.7toJ) O.Olg;^, and ~ —0.3 respectively. The late inter¬ 
nal shock model is disfavored in this burst due to the fact 
that 02 = /? — 2 is signiflcantly smaller than ~ —0.88, un¬ 
less (3 ^ 1.1 which is impossible as indicated by the obser¬ 
vations. Although the derived 02 favors the late external 
shock model, the detected spectral evolution is difficult to 
be explained in this model (see Table 1). 

GRB 050421 — There are two early X-ray flares resid¬ 
ing on a Fh (X t~^'^ light curve. The first stronger flare 
peaks at tpeak = 111 ± 2 s since the trigger, while the 
second weaker flare peaks at t = 154 s (Godet et al. 
2005). The other observed quantities of the first flare are 
Str ~ Stf = 0.07tpeak, and Am,ohs ~ 4 (Godet et al. 2005). 
Because the superposed t~^'^ decaying light curve, which 
lasts at least 500 s since the burst trigger, must have a 
different origin to the two flares, the late external shock 
model for the flare can be ruled out in this burst. For 
the late internal shock model, the intrinsic Am > Am,obs 
requires 0.27 < 2 — (3 < 8.8, which can be easily satis¬ 
fied. Assuming (3 = —1.0, the intrinsic Am of the first 
flare is ~ 50. Therefore the tail emission of this flare at 
late times is about 12 times dimmer than the superposed 
^-3.1 emission, which may originate from the tail of the 
prompt GRB. The value of tb of the first flare can not be 
well constrained. 

GRB 050502B — There are two X-ray flares in this 
burst. Here we just discuss the first very large flare 
peaking at tpeak = 740 s since the trigger. Other ob¬ 
served quantities of this flare are Oi^obs = —<a2,obs = 9.5, 
and Am,ohs ~ 500. The rising and falling times satisfy 
Str I ~ 0.2, Stj^/tpQ^]^ ~ 0.1, and (5t/tpeak I for the 
spike at the peak of the hard band emission (Burrows et 
al. 2005a, b; Falcone et al. 2005). The hardness ratio 
decreases slowly from ~ 1.8 (corresponding to /3 ~ 0.1) 
when t ~ 550 s to ~ 0.7 {(3 ~ —0.4) for t > 850 s and 
has a minor peak at the peak of the hard band (1 — 10 
keV) light curve of the large flare. The measured Am,ohs 
is too large to be explained, since the theoretical maxi¬ 
mum Am,max ~ 33 (when 02 = 3.5) is much smaller than 


this value. The only way to resolve such a problem is that 
the basic assumption on the flare light curve (i.e., eqs. [1] 
& [3]) is oversimplified, as has been directly indicated by 
observations (Ghincarini et al. 2005). The late external 
shock model is ruled out directly in this GRB because 
this basic assumption is quite reasonable, and because the 
pre-flare light curve can be connected with the post-flare 
light curve smoothly by a single power law function of time 
(Fly oc The late internal shock explanation for this 

flare is possible, if the flux density during the whole decay 
phase is described by a pure exponential function of time, 
Fi, = Fjy^pfe exp[—(t — tpeak)/"!"]. The flare increasing fac¬ 
tor in this case is Am = F^_pfc/F^(2tpeak) = exp(tpeak/T). 
Gombining with the observed value, we get the decay time 
scale T « O.lGtpeak = 120 seconds. The equivalent tem¬ 
poral decaying index a2,obs = t/r at t = tpeak and at 
t = 2tpeak is equal to —6.2 and to —12.4, which is consis¬ 
tent with the observed overall decaying index —9.5. The 
ratio Stf/tpeak = (ln2)r/tpeak ~ 0.11 also matches the 
observation. Although the exponential decay assumption 
for the LIS model can explain this large flare, the physical 
reason of this assumption, or equivalently, why the 
tail emission is not detected in this flare, has still not been 
answered. 

GRB 050607 — This burst has been detected with two 
early X-ray flares. The first flare is a weaker one. We 
focus on the second larger flare peaking at tpeak = 310 s 
since the trigger. Other observed quantities of this flare 
are Oi^obs ~ 16, rr2,obs ~ 6.5, St^ Stj^ ~ O.2tp0ak, and 

Am,obs ~ 20 (Pagani et al. 2005). The hardness ratio 
(defined as {H — S)/{H + S) for this burst, Pagani et al. 
2005) decreases slowly from ~ 0.25 (/3 ~ 0.3) at the be¬ 
ginning of the second flare to ~ —0.25 (/? ~ —0.25) at the 
end of this flare. Using a2,obs, Stf/tpeak and Oi^obs, we find 
0,2 ysy —0.6, tb « O.OOtpeak = 29 s, to = 281 s, and ai ~ 1.5. 
The inferred value of 02 favors the late external shock ex¬ 
planation, since in the late internal shock model 02 = i3—2 
is much smaller than —0.6. The inferred is also consis¬ 
tent with the observed temporal index = —0.58±0.07 of 
the post-flare light curve before t ~ 1.2 x 10^ s (Burrows 
et al. 2005b). However, the measured Am,ohs is a little 
larger than the theoretical maximum Am,max ~ 11 when 
02 = —2.4, and even larger than the derived Am = 4.2 us¬ 
ing 02 ~ —0.6 by a factor of 5. This poses a severe crisis 
to the late external shock model explanation. 

4. SEQUENCE OF LATE ACTIVITIES BY CENTRAL 
ENGINES 

The X-ray flares, happened at ~ 10^ s to 10"^ s since 
the trigger of prompt 7-ray emission in the burster’s rest 
frame, are believed to originate from prolonged activities 
of the central engine. Such prolonged activity may be due 
to long-lived intermittent accretion of fragmented materi¬ 
als (or blobs) to the central black hole. Fragmentation may 
take place either in the stellar core of the progenitor during 
its collapse, or in the accretion disk surrounding the black 
hole due to gravitational instability at large radii of the 
disk (King et al. 2005; Perna, Armitage & Zhang 2005). 
Below we discuss possible sequences of central engine ac¬ 
tivities in the context of the LIS/LES models. 

4.1. Gollision between the Early and Late Ejecta 
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Let us first consider that the central engine has been 
active for two periods. This is the basic picture for us to 
understand both prompt GRBs and late X-ray flares. The 
time interval between the beginnings of these two periods 
is denoted as bag. During each activity, the central engine 
ejects one shell whose isotropic-equivalent kinetic energy, 
final mean bulk Lorentz factor and shell width are Ej , Tjo 
and Ajo, respectively. We adopt j = e to denote the early 
ejected shell, while j = I to denote the late ejected shell. 
Each shell may not be uniform and have a variable speed 
ATjo ~ Tjo and variability time scale Atj <C Ajo/c, which 
are required to produce prompt gamma-ray emission by 
internal shocks within the shell. 

In the following, instead of discussing internal shocks 
within the shells, we first study the types of the colli¬ 
sion between the early and late ejecta, which depend on 
whether this collision happens before or after the deceler¬ 
ation of the early shell in the surrounding medium.^ The 
deceleration radius of the early shell is 


Rd 


ec,e 


1.3 X ISM, 

4.0 X 10 ^®ife,53rjQ^2^r^cm, wind. 


(7) 


Subsequently, we define the deceleration time in the 
observer’s frame by tdec.e = 7?dec,e/2rgQC. Gen¬ 
erally, the deceleration time tdec is proportional to 
^i/(3-fc)p-2(4-fe)/(3-fc)^ where k = 0 corresponds to the 
ISM environment while k = 2 corresponds to the wind en¬ 
vironment (n oc R~^). We would like to point out that, 
only the prompt gamma-ray burst originates from internal 
shocks within the early ejecta and the collision between the 
early and late ejecta happens later than the deceleration 
of the early ejecta, and then the arrival time of photons 
emitted from the external shock at i?dec,e driven by the 
early ejecta since the GRB trigger is tdec.e- Usually, there 
are two types of the collision between the early and late 
ejecta. As we show later, if the time interval between these 
two ejecta is large enough (e.g., bag ~ 10® seconds), a new 
type of collision exists, i.e., the late ejecta is decelerated 
when it catches up and sweeps into the non-relativistic tail 
matter of the forward shock driven by the early ejecta. We 
call this type of collision the inner external shock. Below 
we give a quantitative analysis about these three types. 

(i) internal shock — This happens when the late ejecta 
is faster and it collides with the early ejecta before the 
latter is decelerated, i.e., T/g > Ego and i?coi < Rdec,e- 
Here 

Real ~ 2rgQc(bag — Aeo/c) (8) 


by radiation. The property of the subsequent afterglow 
is determined by the property of the merged shell. If 
Ee > El, the merged shell is dominated by the early 
ejecta with Em ~ E^, and T^o ~ Egg. The deceleration 
radius of the merged shell in the circum-burst environ¬ 
ment is Rdec,m ~ Rdec,e, Or equivalently, the deceleration 
time is fdec.m — Rdec,m /2EmgC ~ td ec,e- If Ee < El < 
(E^Q/EgQ)ifg, the energy of the merged shell is Em ~ Ei 
while the Lorentz factor is Emo ~ a/ Ei/E gTeo- In this case 
the deceleration radius and time are Rdec,m ~ Rdec,e and 
tdec.m ^ {Ee/Ei)td ec,e 5 respectively. If (E^p/E^pjAg < Ei, 
the merged shell is dominated by the late ejecta with 
Em ~ El and E^g ~ E;g. The deceleration radius and 
time of the merged shell are Rdec,m ~ Rdec,i > Rdec.e and 
tdec,m ~ tdec,i- From the above analysis, we can see that 
the afterglow may not have the same origin as the prompt 
7 ray emission, if the latter is assumed to originate from 
internal shocks within the early ejecta. 

(ii) refreshed shock — This happens when the time lag 
between the early and late ejecta by the central engine is 
quite large, bag > ^dec,e + Agg/c, or the late ejecta has a 
lower Lorentz factor relative to the early ejecta (Ejg < Egg) 
while the time lag is not large (bag A ^dec.e)- Therefore, 
the collision between these two ejecta happens after the 
deceleration of the early ejecta. The hydrodynamic evo¬ 
lution of the external shock driven by the early ejecta af¬ 
ter its deceleration follows the Blandford-McKee solution 
(Blandford & McKee 1976), 



The external shock, or blast wave, can be approximated as 
a thin shell with a width ARg = i?e/2(3 —b)Eg. When the 
late ejecta catches up with the tail of the thin shell at i?coi, 
the external shock reaches Re = Rcoi + ARe with a Lorentz 
factor Eg « Eg(i?coi)- In the rest frame of the central en¬ 
gine, the time elapsed since the first ejection by the central 
engine is Giapse = Rdec,e/Peoc + ^ dR/fteC, where (3 = 

Vi-i/r^ is the speed in units of c. This time can also be 
estimated by the late ejecta, i.e., Giapse = 7?coi//3/oc + bag- 
Equaling these two elapsed times, we obtain 

7 — 2k 1 1 Rco\ / 3 — k N 

.2(3 - fc)(4 - b)E2(i?ggi) ~ "" V^'^8“4_b^dec.eJ, 

( 10 ) 

or 


is the radius when these two shells begin to collide with 
each other, and ctiag — Agg is the distance from the head 
of the late shell and the tail of the early shell. The second 
condition also reads flag < tdec,e + Agg/c, which requires 
the early ejecta having a very high energy Eg or a rela¬ 
tively low initial Lorentz factor Egg. Because both ejecta 
are cold, this collision is the same as typical internal shocks 
which are responsible for prompt GRBs. 

After the internal shock has finished, the two ejecta 
merge into one shell with energy Em ~ Eg Ei and bulk 
Lorentz factor E„g = a /{Eg + Ei)l(EglY\,^ E Ei/Vf^). 
For simplicity we here neglect the energy loss taken away 


7- 2fc 

3 — k 


tgol 

tdec,e 


r 2 

4 gg 

E 2 
^ (0 



l/(4-fe) 



3- k 

4- fc’ 


( 11 ) 


where Uoi = i?coi/2(4 — A:)Eg(i?coi)c is the observed time 
since the detection of first photons emitted by the early 
ejecta and thus is larger than tdec.e- For Ejg > Egg, since 
the second term on the left side of the above equation can 
always be neglected, the collision happens at 


tgol 


3-k / 

7-2k U'"® 


3 — k ^ / 

-^tdec,e), 


( 12 ) 


^ For simplicity we just consider the thin shell case for the reverse shock propagating into the early shell. The thick shell case could be 
extended in the same way. 
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which is almost independent on the ratio rjo/Feo- This 
means that the time when the two ejecta begin to collide 
detected in the observer’s frame directly measures the time 
lag iiag of these two ejections by the central engine. On 
the other hand, for Tjo <C Feo, the collision happens at 


^col — ^ 

r 7 3-fc F^o 
\7-2k Tl. 

(4-fc)/(3-fe) 

1 ^dec,e 5 ^lag ^lag,criti 

3 — k 

[ 7-2kV'^^ 

3-k \ 

4 _^tdec,ej 

5 ^lag ^ ^lag,criti 

where 

the critical 

value of 

(13) 

flag for a given 

ratio F/o/Fgo is 

^lag,crit — 

[(3 - fc)/(7 - 


iduiu LiQ/LeO i-lag,crit — ~ ~ 

Therefore, if the 


late ejecta moves slower than the early ejecta, we can use 
the collision time tcoi to constrain the contrast between 
the Lorentz factors of these two ejecta when flag is small, 
or to directly obtain flag when it is large. 

The strength of the collision is dependent on the ra¬ 
tio of Fjo to Fe(i?coi)- For the cases of Fjo > Feo and of 
Tio < Feo while flag > bag.crit, we obtain F;o/Fe(i?coi) = 
(F;o/Feo)[(3 - fc)tiag/(7 - 2k)tdec,e]^^~^^ which is 
larger than [(3 — k)/{7 — 2fc)]^/^. On the other hand, 
for the case of F;o < Feo while flag < tiag.crit, we get 
F/o/Fe(i?coi) = [(3 - k)/{7 - 2fc)]^/^. Therefore, a vio¬ 
lent collision requires the initial Lorentz factor of the late 
ejecta much larger than that of the early ejecta, or the 
time interval between these ejecta by the central engine 
long enough (see also Zhang & Meszaros 2002). A for¬ 
ward shock propagating into the downstream fluid of the 
external shock and a reverse shock propagating into the 
late ejecta are developed during this collision. Since the 
total energy of the external shock is increased after the 
two have merged as a whole, the collision is also known 
as the refreshed shock. To produce an obvious variability, 
flattening or bump in the afterglow light curve by the re¬ 
freshed shock requires that the energy in the late ejecta 
must be comparable to or larger than that in the early 
ejecta, or equivalently, in the initial external shock (Ku¬ 
mar & Piran 2000; Zhang & Meszaros 2002). Because the 
very complicated shock jump conditions arises when con¬ 
sidering the forward shock sweeping into a relativistical 
hot medium, the duration of the whole collision must be 
calculated numerically, and it is often found to be much 
longer than tcoi in the observer’s frame. Typical signatures 
from refreshed shocks therefore can not be directly used 
to explain X-ray flares (see figures in Zhang & Meszaros 
2002). Flowever, these signatures are naturally expected 
to appear in afterglow light curves after X-ray flares are 
over and remnants of X-ray flares as late ejecta will collide 
with the ahead external shock. 

(Hi) inner external shock — Above we take the approx¬ 
imation that the gas in the external shock driven by the 
early ejecta is compressed into a thin shell with a single 
bulk Lorentz factor. In fact, according to Blandford & 
McKee (1976), the distributions of the Lorentz factor 7 bm 
and the number density ubm (in the observer’s frame) of 
shocked media can be described as simple functions of the 
self-similar variable y, 

7BM = 2-i/2f,x-i/2, ubm = 2F2x-(7-2'=)/(4-'=)„, (14) 

where x = 1 -I- 2(4 — A:)(Ai?/i?e)Fg, Ai? = — re is the 


radial distance of the shocked media to the shock front, 
and n = AR~^ is the circum-burst media density. Far be¬ 
hind the shock front, say rg < i?e/ 2 , a uniform cavity or 
bubble is shaped with the number density, 


^bub 


3.7 X ISM, 

5.6 X wind. 


(15) 


Note that this bubble is cool and quasi-static. The shocked 
media with 7 bm — 1 is located at re,„ = [(15 — 4A:)/(16 — 
4fc)]i?e- The total mass in the non-relativistic tail of the 
external shock is 


Me,n = f ^irr'^UBumpdr = CfcF2('=-3)/(4-fc)Mgw, 

(16) 

where Mgw = 47 rnmpi?g /(3 — k) is the mass of swept cir- 
cum burst media, and the coefficient Ck = 0.96 for /c = 0 
and Ck = 0.43 for k = 2. 

Suppose that the late ejecta reaches rg_„ while still not 
being decelerated in the non-relativistic tail of the exter¬ 
nal shock. In the rest frame of the central engine, the 
time elapsed since the first ejection by the central en¬ 
gine is teiapse — 77g[gg g//?gQC 4- dRj(deC. It is also 

equal to <eiap se — ^e,n lldioc + flag. Assuming F/o, Fg 
much larger than unity, we have rg^„ ~ (15 — 4fc)c[tiag — 
(3 — fc)tdec,e/(4 — k)], or the external shock radius Re — 
4(4 - A:)c[tiag - (3 - A:)tdec,e/(4 - k)]. Since Re > Rdec,e, 
the time lag between the two ejections is required to be 
very long, 

tiag >Fydec,e/[2(4-A:)]. (17) 

For example, adopting Fgo ~ 100 and tdec.e ~ 100 s, we 
require flag > 1 day. This long lag challenges the most 
popular central engine models for GRBs, however, it is 
supported by recent observations (e.g., Fox et al. 2005 
for an X-ray flare happened at 16 days after the short 
burst GRB 050709). Now we turn to derive the critical 
value of F/o, i.e., F/o,crit- A late ejecta with this criti¬ 
cal F/o will be decelerated exactly at Rdec,i = where 
i?dec,i is the deceleration radius of the late ejecta. Using 
^i/rfo.critC^ = ^e,n and Ug « MswF^c^, we get 


F/o,c 


VUg/ ® 

tdec,e 

\Ee) L2(4-fc)tlag- 


(3-fc)(7-2fc)/[2(4-fc)] 

(18) 


where Ck is ignored for simplicity. Late ejecta with same 
flag but F/o > F/o,crit will be decelerated before reach¬ 
ing re,n, i-e., Rdec,i < 're,n- Therefore, the external shock 
driven by the late ejecta and propagating into the non- 
relativistic tail of the outer external shock are called the 
inner external shock. Since the mass in the bubble within 
Rdec,i is equal to Mg_„(i?dec,i/?'g,n)^, the deceleration ra¬ 
dius of the late ejecta is i?dec,i = ''’e,n(F^i/Mg^„FfpC^)^/^. 
Again, letting the elapsed time Rdec e(PeocE fp' dR/PeC 

equal to Rdec,i /Pioc+tug, we obtain the radius of the outer 
external shock when the late ejecta has been just deceler- 
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ated, 


Re 




2(4-fc)/[(3-fc)(7-2fc)] 


if Tio ,crit < Tio < DkTiQ ,crit 5 




if -DfcFio jCrit < Fjo, 


(19) 


where Dk = (15 - 4fc)3/2(16 - 4fc)(2fc"-iofc+9)/(8-2fc) jg 
~ 1.3 X 10^ for fc = 0 and ~ 3.9 for k = 2. The dis¬ 
tance between the outer external shock and the late ejecta 
when the latter reaches its deceleration radius is always 
Re — Rdec,i — ctiag. Similar to the outer external shock, 
the X-ray emission from the inner external shock peaks at 
the moment of its deceleration. The time scale of the peak 
in the observer’s frame is ~ tdec,i = 7?dec,i/2F^QC, which is 
much shorter than the lag time <iag between the inner and 
outer shocks. Therefore, the observed emission from the 
inner external shock behaves as a very short spike. The 
inner external shock will finally catches up with the main 
part of the outer external shock. Because the gas just be¬ 
hind each external shock are both hot, strong forward and 
reverse shocks are hardly to be developed when they begin 
to collide. 

In Fig. 3 we visualize the above three types of collision 
between the early and late ejecta. Assuming the properties 
of the early ejecta are known, the type of collision depends 
mainly on the ratio of the initial Lorentz factor of the late 
ejecta to that of the early ejecta, F/o/Fgo, and on the time 
lag flag between these two ejections by the central engine. 
More details can be found in the caption of Fig. 3. Here 
we note that for the refreshed shock, the relative speed of 
the late ejecta in the rest frame of the outer external shock 
must be supersonic, i.e., larger than This condition 

is equivalent to that the relative Lorentz factor, which is 
equal to [Fio/Fe(i?coi)+ Fe(i?coi)/F/o]/2 as long as F/o and 
Fe(7?coi) much larger than unity, must be greater than 
1.22 (Zhang & Meszaros 2002). Therefore the criterion 
for forming a refreshed shock is F;o/Fe(i?coi) > 1-92. We 
can see from Fig. 3 that if the refreshed shock is strong, 
the observed beginning time of the collision directly scales 
to the time lag by tcoi — [(3 — A:)/(7 — 2fc)]tiag. This is 
especially important because only strong refreshed shocks 
can be easily detected by a brightening or flaring signature 
in the light curve. In the following we just consider the 
strong collision cases, including the internal shock case 
(F/o > Feo and flag < tdec,e)j and the strong refreshed 
shock cases (F/o > Feo and t/ag > tdec,e, or F/o < Feo and 

^lag ^ flag.crit)- 


4.2. Types of X-ray light curves 

Now let us turn to discuss the entire X-ray light curve 
from the early and late ejecta. Temporal variabilities in 
at least one of the two ejecta are required to produce the 
prompt GRBs. There are four basic types of X-ray light 
curves, depending on the prompt GRB emission originat¬ 
ing from the early ejecta or from the late ejecta, as well 
as the collision between the two ejecta taking place earlier 
than the deceleration of the early ejecta or not. 

Type 1 — Prompt gamma-ray burst arises from inter¬ 
nal shocks in the early ejecta, and the collision between 
the early and late ejecta happens earlier than the decel¬ 
eration of the early ejecta, i.e., i?coi < 7?dec,e- This case 


corresponds to F/o > Feo, and t/ag < tdec,e- 

The prompt GRB is produced at the internal shock 
radius of the early ejecta, i?int,e ~ 2FgocAte- If there 
is also time variability in the late ejecta, then a second 
burst will be produced by the internal shock within the 
late ejecta at Rint,i ~ 2 F^QcAt/. At this moment we have 
the equation for the elapsed times in the burster’s frame, 
Rint,e//?eOC “f (-^GRB R'wit^e )/^— ^lag 4“ Rint,l/Pl0^j where 
7?grb is the radius of prompt gamma-ray photons when 
the internal shock of the late ejecta happens. The time 
interval between the beginnings of the two bursts, or the 
to of the second burst, is 

toil) = ~ = ^lag + At/ - Ate - tlag. ( 20 ) 

When Rint,i > Rint,e, the main emission of the second 
burst is in the X-ray band and the internal shock of the late 
ejecta produces the first X-ray flare, providing Af/ <C t/ag- 
Since the light curves of most of observed X-ray flares are 
less structured, the late ejecta is required to be composed 
of a few (two or three) mini-shells. A second X-ray flare is 
unavoidable because of the collision between the early and 
late ejecta. The equation for the elapsed times at this col¬ 
lision is Rint,e//4e0^ 4“ (RgRB 7?int,e)/c — t/ag 4- Rcolf 
and i?GRB here is the radius of prompt gamma-ray pho¬ 
tons when the two ejecta collide with each other. The to of 
the second X-ray flare relative to the trigger of the prompt 
GRB is 

toic) = « (l + ^)ilag - ^TgrB, ( 21 ) 

^ ^ lo' ^10 

where Tgrb ~ Aeo/c is the duration of GRB. After the 
collision, the merged shell with energy Em and Lorentz 
factor Fmo sweeps into the circum-burst medium and will 
be decelerated at Rdec,m- The maximal flux of X-ray emis¬ 
sion of the afterglow is also produced at this radius. By 
using the equation for the elapsed times at the decelera¬ 
tion, R'lnt^e/Ped^ 4- (7?GRB R'mt^e)/C- — ^lag 4- Rco\/Pw^ 4- 
(7?dec,m — Rco\)/Pmoc, we have the peak time of the X-ray 
afterglow in the observer’s frame. 


ipk(AG) 


-^GRB Rdec^m 
C 



( 22 ) 


where we have assumed Aeo/c <C tiag- As discussed in the 
last subsection, the property of the merged shell is deter¬ 
mined by the energy ratio Ei/Ef.. If Ei < Eg, since Fmo ~ 

FgO and tdec^m ^ ^dec,e ^ ^lag; ll&VG tp]^(AG) ~ t(\ ec,m 5 

which means that the to effect can be neglected in this 
case. The prompt GRB and afterglow have the same ori¬ 
gin. If Ee < El < (Ffp/F2p)F;g, since r rnO ^ 
and ^dec,m (-^e/^z)^dec,e; WG liaVG tpk(AG) ~ (1 
Ee/Ei)t\a,g 4- tdec,m- The X-ray afterglow would behave 
as an X-ray flare only if tdec,m 'C (1 — E^/Ept^g, which 
requires EijE^ > tdec,e/tiag and F/o/Feo > i/tdec.e/tiag- 
Inversely, the X-ray afterglow does not suffer from the 
to effect. If {rfg/TlQ)Ee < Ei, since F^o ~ F/o and 
tdec,m ^ tdec.h have tpk(AG) ~ ^lag 4- tdec,m- The 
criterion for the X-ray afterglow to be an X-ray flare is 



tdec.m < ilag, Or Ej/Ee < (ilagAdec,e'"(r/o/reo)® 

and Fjo/reo ^ \/idec,e/^iag- In conclusion, the late ex¬ 
ternal shock model for X-ray flares requires the energy 
of the late ejecta is much larger than the energy of the 
early ejecta, and the ratio of initial Lorentz factors sat¬ 
isfies F/o/Feo ^ a/ tdec,e/tiag- This model also implic¬ 
itly requires a very low initial value of Feo, because of 
tdec.e > tlag ~ 10® — 10"^ SeCOnds. 

Type 2 — Prompt gamma-ray burst arises from internal 
shocks in the early ejecta, and the collision between the 
early and late ejecta happens later than the deceleration 
of the early ejecta, i.e., i?coi > ^dec,e- We just consider 
a strong refreshed shock being developed during the col¬ 
lision. This corresponds to the case of F/o > Feo, tiag > 

^dec,ej Cr the case of F/o ^ F^g? ^lag ^ ^lag,crit ^ ^dec.e- 

Similar to the analysis in Type 1, the prompt gamma- 
ray emission is produced by internal shocks in the early 
ejecta at i?int,e ~ 2FggcAte- In the observer’s frame, the 
X-ray afterglow from the external shock driven by the early 
ejecta will peak at 

ipk(AG) = tdec,e- (23) 

If a second burst is produced by the internal shock within 
the late ejecta, then the to for this burst is 

to{l) « flag, (24) 

which means the second burst happens later than the 
deceleration of the early ejecta in the observer’s frame. 
Subsequently, when the late ejecta begins to collide with 
the external shock originating from the early ejecta, the 
equation for the elapsed times is (i?GRB — -Rint,e)/c = 

(.I^dec.e .I^int,e)//?eOC“t “^ dR/ /?eC, where R^ — i?col[l“t“ 

1/2(3 —A:)Fg(i?coi)]- The to of the refreshed shock emission 
relative to the prompt gamma-ray emission is 

. 7 — 2k 3 — k , , 

fo(c) ~ l^tcQi ^tdec,e — ^lag- ( 2 o) 

This means that the refreshed shock emission overlaps the 
internal shock emission from the late ejecta in the ob¬ 
server’s frame. 

Type 3 — Prompt gamma-ray burst arises from inter¬ 
nal shocks in the late ejecta, and the collision between the 
early and late ejecta happens earlier than the deceleration 
of the early ejecta, i.e., i?coi < -Rdec.e- This corresponds 
to the case of F/o > Feo, tiag < ^dec.e- Note that in this 
case the internal shock in the early ejecta is forbidden, 
otherwise this emission would precede the emission from 
internal shocks in the late ejecta by t/ag in the observer’s 
frame and would trigger the detector. In other words, 
there should be no time variability in the early ejecta. 

After producing the prompt GRB at ~ 2F/pcAf/, 

the late ejecta will collide with the early ejecta at the ra¬ 
dius Rco\- The equation for the elapsed times at the colli¬ 
sion is i?i„t.///3/oc-l- (i?GRB - Rint,i)/c = Real/Pioc. There¬ 
fore, the beginning time to of this collision is 

fo(c)« S^(<lag-Aeo/c), (26) 

^ 10 

which is smaller than the lag tug. When the collision is 
over, the early and late ejecta merge into a shell with en¬ 
ergy E^ and Lorentz factor F^o- The emission of X-ray 


afterglow peaks when the merged shell reaches to its decel¬ 
eration radius Rdec,m- The equation for the elapsed time 
at the deceleration is Rint,i/^ioc + {Rgrb - Rint,i)/c = 
Rcoi/Pioc + (i?dec.m “ Rcoi)/f^moC- The peak time of the 
X-ray afterglow in the observer’s frame is 

tpk(AG) « tdec.m - - 3^) (ilag " ^eo/c). (27) 

Below we neglect the Ago /c term for simplicity. If 

El < since BgO S-nd tdec,m ^dec,e ^ ^lag 5 

we have ^pk(AG) ~ ^dec,m ^lag ^dec,m • If Ee < 

El < (Y^q^E^^ since ~ \/ EijE^Y eo and tdec,m 

(_£/e/-£'/)^dec,e 5 We have ipk(AG) ~ (B/e/A'/)(tdec,e ^lag) 

^dec,m ■ If (F/o/F^o)^e < since F^o ~ F/o and 

tdec,m ~ tdec,i, we have tpk(AG) « tdec.m- In conclusion, 
the to effect for the light curve of X-ray afterglow can be 
neglected in this case. 

Type 4 — Prompt gamma-ray burst arises from inter¬ 
nal shocks in the late ejecta, and the collision between the 
early and late ejecta happens later than the deceleration 
of the early ejecta, i.e., i?coi > Rdec,e- We just consider 
a strong refreshed shock being developed during the col¬ 
lision. This corresponds to the case of F/o > Feo, tiag > 
tdec,e^ Cr the case of F/o ^ FgO, t/ag ^ tlag,crit ^ f'dec,e' As 
in Type 3, the internal shock in the early ejecta is forbid¬ 
den. 

The equation for the elapsed time when the early ejecta 
reaches its deceleration radius is Rdec,e/Peoc = tiag + 
Rint,i/f3ioc + {Rgrb - Rint,i)/c. The X-ray afterglow of 
the external shock driven by the early ejecta peaks at 

^pk(AG) = tdec,e ^lag At/ ~ tdec,e ^lag- (28) 

This means that the peak emission of X-ray afterglow hap¬ 
pens before the main gamma-ray burst in the observer’s 
frame. Contrary to the above three types, the beginning 
time of this kind of X-ray afterglow relative to the prompt 
burst is negative, to = —tiag- The detected light curve of 
the X-ray afterglow behaves as 



where t is the time in the observer’s frame since the trig¬ 
ger of the prompt burst, and a 2 ~ ~1 is the intrinsic 
temporal decaying index. The observed temporal index 
is Q; 2 ,obs = [t/{t + tiag)]c( 2 , which approaches zero when 
t tiag while approaches a 2 when t > flag. Such a nega¬ 
tive to could explain the slowly decaying early X-ray after¬ 
glows recently discovered by the Swift satellite (Zhang et 
al. 2005; Nousek et al. 2005). Subsequently, when the late 
ejecta begins to collide with the external shock originating 
from the early ejecta, the equation for the elapsed times 

is t/ag T Rinty/PlO^ T (Rgrb Rint,/)/^ — Rdec,e//^eOC “t“ 

iRdec e where Re = Rcoi[l + 1/2(3 - A:)Fg(Rcoi)]. 

The to of the refreshed shock emission relative to the 
prompt gamma-ray emission is 

.. 7 — 2k 3 — k , , 

fo(c) ~ ^tcol ^lag “t” ~jTdec^e ~ 0. (30) 
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This means that the refreshed shock emission happens at 
nearly the same time with the prompt GRB produced by 
the internal shocks within the late ejecta in the observer’s 
frame. 

Fig. 4 shows the above four types of X-ray light curves. 
Some caveats must be given here. First, for Types 1 Sz 2, 
the internal shock within the late ejecta (thin solid lines) 
may not take place if there is no temporal variability in 
the late ejecta. Secondly, refreshed shocks are assumed to 
be strong in the hgure {Types 2 k, 4)- Weak or mildly 
strong refreshed shocks will hardly influence the X-ray 
light curve. Thirdly, the temporal behavior of afterglow 
emission in Types 1 light curve depends on the property of 
the merger of the two ejecta. It may behave as an X-ray 
flare if the property of the merger is determined by the 
late ejecta with Ei > and Fjo/Feo > (tdec,e/iiag)^/^- 
This corresponds to the late external shock model. Other¬ 
wise the afterglow light curve is normal and its beginning 
happens when the prompt GRB is produced. 

Now we begin to classify the observed X-ray light curves 
by using the above results. As discussed in Section 3, GRB 
050406 has one X-ray flare with tpeak = 213 s. The de¬ 
rived intrinsic is ~ —0.56. We note that the post-flare 
light curve behaves as Ei, oc i-i-58±o.i7 £qj. ^ ^ g 
El, oc f:-o-50±o.i4 ^ ^ ^ g (B^-pows et al. 2005; 

Romano et al. 2005). In principle, there are two possible 
explanations for this burst. In the first explanation, the 
X-ray light curve of GRB 050406 belongs to the Type 1 
light curve. The X-ray flare is caused by the late external 
shock. The intrinsic is consistent with the temporal 
index of the very late afterglow light curve. In the second 
explanation, the X-ray light curve belongs to the Type 2 
light curve. The X-ray flare is interpreted as the internal 
shock emission within the late ejecta. The deceleration 
time of the early ejecta satisfies fdec.e < 100 s and the af¬ 
terglow decays as E^ ex The collision between the 

late ejecta and the external shock must be weak or mildly 
strong because there is no obvious signature in the light 
curve just after the X-ray flare. In this explanation, we 
also need a continuous energy injection from t ~ 4 x 10^ s 
to explain the flattening of the light curve. 

GRB 050421 definitely belongs to the Type 1 light curve, 
because the two X-ray flares reside on a Fiy oc emis¬ 
sion background, which is the tail emission of the prompt 
gamma-ray burst. The first X-ray flare is naturally inter¬ 
preted as from the internal shock within the late ejecta, 
while the second X-ray flare is from the subsequent in¬ 
ternal shock between the late ejecta and the early ejecta. 
It is intriguing of this GRB that the tail emission lasts 
for about 10^ s and the afterglow emission has not been 
detected. 

The X-ray light curve of GRB 050502B is classified 
to be the Type 2 light curve. The afterglow decays as 
El, oc ® with the deceleration time tdec.e < 50 s. The 
Hrst X-ray flare peaks at fpeak ~ 740 s. The fluence of 
this X-ray flare is comparable with that of the prompt 
GRB. This means the energies contained in the late ejecta 
and in the early ejecta are comparable. One may expect 
a strong refreshed shock signature nearly simultaneously 
following this X-ray flare would be detected. However, 
a second X-ray flare has been detected, but it happened 
at t ~ 4 X 10"^ s, much later than the first X-ray flare. 
This implies the refreshed shock is weak, which requires 


the initial Lorentz factor of the late ejecta is smaller than 
that of the early ejecta, F;o < Feo, and the time lag 
is no longer than the critical value, bag < bag.crit, or 
Ro/reo < Taken tdec,e ^ 30 S 

and bag ~ 300 s, we have Fjo/Feo < 0.42 for b = 0 (ISM) 
and F/o/Feo < 0.56 for b = 2 (wind). The second X-ray 
flare happens at f ~ 0.5 day requires another late ejecta 
by the central engine to catch up with the external shock 
and produce a strong refreshed shock. 

GRB 050607 may belong to the Type 1 light curve. The 
Hrst weaker X-ray flare can be either from the internal 
shock within the late ejecta, or from the internal shock 
between the early and late ejecta. The second large X- 
ray flare is interpreted as the late external shock emission. 
However, as discussed in Section 3, the main disadvantage 
of this late external shock explanation is the theoretical 
flare increasing factor Am. ~ 4 is much smaller than the 
observed Am,ohs ~ 20. 

5. DISCUSSION AND CONCLUSIONS 

In this paper we have quantitatively analyzed late X-ray 
flares in the frameworks of the late internal shock model 
and late external shock model. As we have shown above, 
both the late internal shock model and late external shock 
model require late time activities of central engines. Some 
of previous works suggested the late internal shock origin 
for X-ray flares (Burrows et al. 2005a; Zhang et al. 2005; 
Fan & Wei 2005), while others suggested the late external 
shock origin (Piro et al. 2005; Galli & Piro 2005). In fact, 
these two kinds of late shocks can coexist within a certain 
gamma-ray burst. GRB 050607 may be such a possible 
candidate. 

Here one caveat should be made for the late external 
shock model. An X-ray flare originating from the late ex¬ 
ternal shock comes to being only if the prompt gamma-ray 
burst is produced by internal shocks within the early time 
ejecta, and the late time ejecta carries more energy than 
the early one and its initial Lorentz factor F;o must be 
larger than the initial Lorentz factor of the early ejecta 
Feo by a factor of (tdec,e/bag)^^^ > 1 , where bag is the 
time separation between the two ejecta by the central en¬ 
gine and tdec.e is the deceleration time of the early ejecta 
in the circum-burst medium. The merger of the two ejecta 
is mainly dominated by the late ejecta with a much shorter 
deceleration time relative to tdec,e- Therefore the begin¬ 
ning of the afterglow from the merger happens much later 
than the prompt GRB and the peak of afterglow behaves 
as an X-ray flare. Note this description of the late inter¬ 
nal shock model is different from that in Piro et al. (2005) 
and Galli & Piro (2005). They attributed the late external 
shock to arise from a thick shell which is ejected by the 
long active central engine. Despite different descriptions, 
this model requires the to effect in the afterglow emission. 
Given to ~ bag) the smaller the deceleration time of the 
merger tdec.m) the more remarkable the afterglow as an X- 

ray flare. Since tdec.m oc a large 

energy or especially a large Lorentz factor will easily cause 
an X-ray flare. This result has also been obtained by Galli 
& Piro (2005). 

Theoretically, there are four basic types of X-ray light 
curves if the central engine has two periods of activities 
before it entirely ceases. According to our analysis, the 
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observed X-ray light curves of GRBs with X-ray flares 
all belong to Type 1 or Type 2. This means the prompt 
gamma-ray emissions of these GRBs result from the early 
ejecta, while the late X-ray flares are produced by the in¬ 
ternal shock within the late ejecta, by the internal shock 
between the early and late ejecta, or by the late exter¬ 
nal shock. The difference between Type 1 and Type 2 is 
that for the former the collision between the early and late 
ejecta happens before the deceleration of the early ejecta, 
while for the latter the situation is inverse. One may won¬ 
der where are Type 3 and Type 4 X-ray light curves? In 
these two types, the prompt gamma-ray burst is produced 
from the late ejecta. There would be no internal shocks 
within the early ejecta, or these internal shock emission is 
too weak to be detected. It is also possible that the inter¬ 
nal shock emission within the early ejecta is very weak and 
regarded as the precursor of the main burst. According to 
our classification. Type 3 X-ray light curves correspond to 
the case that the collision between the early and late ejecta 
takes place before the deceleration of the early ejecta. The 
emission from the collision has a small time lag relative to 
the prompt gamma-ray emission and can be regarded as 
the last pulse of the prompt GRB. The afterglow emis¬ 
sion in Type 3 hardly suffers the to effect. Therefore it 
is hard to distinguish Type 3 X-ray light curves from the 
light curves in which the central engine has only one ac¬ 
tive period. However, Type 4 X-ray light curves may have 
already been detected. Some of recently discovered X-ray 
afterglows whose initial decay is extraordinarily slow may 
originate from this type. 

In this paper we also investigate a new kind of colli¬ 
sion between the early and late ejecta, the inner external 
shock. Such a shock will be developed when the late ejecta 
is decelerated by the non-relativistic tail of the outer ex¬ 
ternal shock, which is driven by the early ejecta. To de¬ 
velop the inner external shock requires the time lag flag 
between these two ejecta is very long, typically > 1 day, 
and the baryon loading of the late ejecta is very low, or 
equivalently, its initial Lorentz factor must be larger than 
a critical value. Since the required bag is not very ex¬ 
treme compared to recent discoveries and the low baryon 
loading is also plausible, the inner external shock maybe 
exist. Because the time scale of emission from this shock 
in the observer’s frame is equal to its deceleration time and 
therefore is much shorter than flag, the emission behaves 
as a short spike. Detecting such a spike in the typical 
afterglow time scale {t ~ bag ~ 1 day) is quite difficult. 
We therefore just make a prediction of the inner external 
shock emission in this paper, and have not considered such 
kind of emission in the above four basic X-ray light curves. 
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For simplicity we have only considered the central en¬ 
gine having two periods of activities. In reality the central 
engine may have more than two periods of activities. The 
multiple X-ray flares detected in the high redshift gamma 
ray burst GRB 050904 indicate that the central engine of 
this burst has been active for many periods (Watson et al. 
2005; Gusumano et al. 2005). Cusumano et al. (2005) 
reported that the hardness ratio for t > 100 s is nearly a 
constant, iJ(14 —73keV)/S'(1.4—14keV) ~ 0.2, which cor¬ 
responds to ~ —0.9. The X-ray light curves from mul¬ 
tiple ejections by the central engine are complicated and 
difficult to be classified. As for GRB 050904, Zou, Xu & 
Dai (2005) found that ordered late internal shocks formed 
by collisions between late ejecta and the earliest ejecta 
can reproduce the temporal evolution of peak luminosity 
of X-ray flares and the entire X-ray spectral evolution. 

At last we would like to point out other ways besides 
those developed in this paper (see Section 3) to distinguish 
an X-ray flare caused by a late external shock from other 
X-ray flares in the same GRB by late internal shock(s). As 
we know in most pulses of GRBs it is found that higher 
energy photons arrive to the observer earlier than lower en¬ 
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lag phenomenon. Fenimore et al. (1995) found that the 
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at higher photon energy E, i.e., W oc This is called 

the spectral-width relationship. If these two phenomena 
have been detected in a particular X-ray flare, then the 
X-ray flare can be proved to originate from a late inter¬ 
nal shock. Statistics help us to understand X-ray flares. 
In prompt GRBs, several relationships have been found, 
such as the isotropic gamma-ray released energy Eiso and 
spectral peak photon energy Ep relationship (Amati et al. 
2002 ), the peak luminosity Lp and spectral peak photon 
energy Ep relationship (Yonetoku et al. 2004), the peak 
luminosity Lp and spectral time lag r relationship (Nor¬ 
ris, Marani & Bonnell 2000). Moreover, Liang, Dai & 
Wu (2004) found there exists a relationship between the 
flux and Ep within a GRB with time resolved spectra (see 
also Borgonovo & Ryde 2001). Therefore, investigating 
the above relations in X-ray flares is urgent since it can 
provide direct evidence that whether late X-ray flares and 
prompt GRBs have the same mechanism. 
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Table 1 

Temporal and spectral indices (F^ oc of the forward shock emission in ism 


regime 

a 

t < tb (NRS) 

ai 

t < tb (RRS) 

CHi 

[t > tb) 

0(2 

Sij'p j f 

(NRS vs RRS) 

V <Vc<Vm 

1/3 

11/3 

4/3 

1/6 

none 

V <Vm<l'c 

1/3 

3 

4/3 

1/2 

none 

Vc < V < Vm 

-1/2 

2 

1/2 

-1/4 

0.02 vs 0.05 

< ly <iyc 

-(p-l)/2 

3 

(3-p)/2 

-3(p-l)/4 

0.18 vs 0.71 

max{iyjn, i^c} < v 

-p/2 

2 

cT 

1 

3 

T 

-(3p-2)/4 

0.35 vs none 


Note.—T he indices are derived under the assumptions that the hydrodynamics is adiabatic and 
the emission is from the synchrotron radiation of shock-accelerated electrons, whose initial energy 
distribution follows dN/dje oc The ratio 5tr/Stf is estimated by assuming the observing band 
in the same spectral regime around tf, and adopting p = 2.2. 
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Table 2 

Temporal and spectral indices cx of the forward shock emission in stellar wind 


regime 

p 

Qfl {t < tb) 

a2 (t > tb) 

Stj' j f 

Vc<V<Vm 

-1/2 

1/2 

-1/4 

0.05 

Vm <V <Vc 

-{p-l)/2 

-(p-l)/2 

-(3p-l)/4 

none 

max{i/m, Vc} < V 

-p/2 

cT 

1 

1 

-(3p-2)/4 

none 


Note. —Similar to Table 1 except for a stellar wind environment. The 
theoretical ai in the case of is also presented in the table, 

although Vc is usually below Vm and the X-ray frequency for typical 
parameter values. 



Fig. 1.— The light curve in grey describes the intrinsic temporal evolution of the emissions from one internal or external 
shock. The light curve in black peaked at tpeak is the time-shifted intrinsic light curve by to- The parameters are defined 
in the text. 
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photon index r=p-1 


Fig. 2.— Hardness ratio (ratio of the hard-band count rate to soft-band count rate) as a function of the effective photon 
index F, assuming the 0.2 — 10 keV spectrum can be fitted by one single power-law Fi, (x . 
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ISM wind 



Fig. 3.— Three types of collision when the late ejecta catches up with the early ejecta, depending on the ratio of T/o to 
Teo and tiag (which is scaled to the deceleration time of the early ejecta). Left and right panels correspond to two kinds 
of circum burst environment, i.e., the interstellar medium (left) and the free wind (right). The upper left grey region in 
each figure denotes the region for internal shocks. We simply neglect the width Aeo of the early ejecta. The right grey 
region (plotted assuming Teo = 100, Ei = O.lEe) denotes the inner external shock region, in which the late ejecta can be 
decelerated in the non-relativistic tail of the early external shock. The rest region is for refreshed shocks. Different solid 
lines correspond to different values of tcoi/tdec,e, ranging from 1 to ^ 10^ as marked beside each line. Different dotted 
lines correspond to different values of r/o/re(i?coi), which are marked on the lines. If r/o/re(i?coi) < 1.92, the refreshed 
shock is suppressed by the thermal energy in the outer external shock. 
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Fig. 4.— Four types of X-ray light curves, depending on the prompt GRB produced by internal shocks within the early 
ejecta {Types 1 h 2) or by internals shocks within the late ejecta {Types 3 Sz 4)- In each case, the X-ray light curve 
is further determined by whether the collision between the two ejecta happens earlier than the deceleration of the early 
ejecta {Types 1 & 5) or not {Types 2 Sz 4)- Thick solid lines represent the X-ray emission during the prompt GRB phase. 
The decaying index of the tail emission of the prompt GRB is—2-|-/9~—3. Thin solid lines correspond to the X-ray 
emission produced by the internal shock of the late ejecta. Dashed lines correspond to the emission from the collision 
between the two ejecta, which can be either the direct internal shock (Types 1 Sz 3) or the refreshed shock (Types 2 Sz 4)- 
Dotted lines correspond to the emission from the external shock. 








